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Cosmic Ray composition around the knee from EAS electromagnetic and 
muon data 

The EAS-TOP and MACRO Collaboration * 


We report the analysis of the Ng-N^ coincident events collected by the MACRO/EAS-TOP collabo¬ 
ration at the Gran Sasso Laboratories. The result points to a primary composition becoming heavier 
around the knee of the primary spectrum (in the energy region 10^® — 10^® eV). The result is in very good 
agreement with the measurements of EAS-TOP alone at the surface, wich detects muons with energy 
E^> 1 GeV and uses the same (QGSJET) interaction model. 

1. Introduction and the detectors 

The study of the primary composition in the Extensive Air Shower energy region requires the use 
of different observables in order to cross check the information and reduce the dependence on the in¬ 
teraction model and propagation codes used. At the National Gran Sasso Laboratories a program of 
exploiting the surface shower size measurements from EAS-TOP (2005 m a.s.l.), and the high energy 
muon measurements = 1.3 TeV) performed in the deep underground laboratories (MACRO) has 
been developed. Such muons in fact originate from the decays of mesons produced in the hrst interactions 
in the atmosphere and from a quite different rapidity region than the GeV muons usually used for such 
analysis (xf > 0.1 or 0.2). 

The two experiments operated in coincidence for a live time of AT = 960.12 days between 1990 and 
2000. The number of coincident events collected with the two detectors operating in the final configuration 
amounts to 28160, of which 3752 have shower cores inside the edges of the array (’’internal events”) and 
shower size Ne > 2.10^, and 409 have Ne > 10® ®^, i.e. above the observed knee at the corresponding 
zenith angle. We present here an analysis of the full data set in terms of the QGSJET interaction model 
as implemented in GORSIKA. For a comparison of different interaction models see |0]. 

The EAS-TOP array is located at Gampo Imperatore (2005 m a.s.l., « 30° with respect to the vertical 
of the underground Gran Sasso Laboratories, corresponding to 930 gr cm“^ atmospheric depth). Its e.m. 
detector (to which we are mainly interested in the present analysis) is built of 35 scintillator modules 
10 m^ each, including an area A « 10® m^. The array is fully efficient for Ng > 10®. Its reconstruction 
capabilities of the EAS parameters (for internal events) are: ^ 10% above Ng w 10®, and A9 ft! 0.9° 

for the EAS arrival direction. The array and the reconstruction procedures are fully described in 

MAGRO, in the underground Gran Sasso Laboratory at 963 m a.s.l., with 3100 m w.e. of minimum 
rock overburden, is a large area multi-purpose apparatus designed to detect penetrating cosmic radiation. 
The lower part of the MAGRO detector has dimensions 76.6 x 12 x 4.8 m®. A detailed description of the 
apparatus can be found in ||]. 

In this work we consider only muon tracks which have at least 4 aligned hits in both views of the hori¬ 
zontal streamer tube planes over the 10 layers composing the whole detector. The standard reconstruction 
procedure of MAGRO |Q] has been used. 

The two experiments are separated by a thickness of rock ranging from 1100 up to 1300 m, depending 
on the angle. 

The muon energy threshold at the surface for muons reaching the MAGRO depth ranges from E^J^ = 1.3 
TeV to Ejl^ = 1.8 TeV inside the effective area of EAS-TOP. Event coincidence is established off-line, 
using the absolute time given by a GPS system with an accuracy better than 1 fis. 

Independent analyses of the two arrays are in ||] and . 

*The complete list of authors will be reported in the summary of the TAUP2001 Proceedings volume. 
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2. Analysis and results 

The analysis technique has to be adapted to the specific trigger requirements (both surface and under¬ 
ground detectors fired) with defined acceptance area for the EAS array (internal events), but undefined 
for the underground one. Therefore, the main experimental feature to be considered is the muon mul¬ 
tiplicity distribution in different intervals of shower sizes. We have chosen six intervals of shower sizes 
around the knee: 

5.20 < Logw{Ne) < 5.31; 5.31 < Logio{Ne) < 5.61; 5.61 < Logio{Ne) < 5.92; 5.92 < Logw{Ne) < 6.15; 
6.15 < LogioiNe) < 6.35 and 6.35 < LogiQ^Ne) < 6.7. 

Within each size bin the muon multiplicity distribution has been fitted with a superposition of pure p 
and Fe components, or light (L) and heavy (i7) admixtures containing equal fractions of p and Fie or 
Mg and Fe respectively. All spectra in the simulation have slope 7 = 2.62. 

The fit has been performed in the quoted six windows by minimizing the following expression for each 
multiplicity distribution: 

^,2 ^ (ivr-Piy-P2ivn^ 

j ^i,exp 

where in the number of observed events in the i-th bin of multiplicity, Np {N^) and are 

the number of simulated events in the same i-th multiplicity bin from the p {L) and Fe (H) components, 
respectively; pi and p 2 are the parameters (to be fitted) defining the fraction of each mass component 
contributing to the same multiplicity bin. The results of the fits have been normalized to reproduce the 
observed number of coincident events in each size bin. 

To infer the mass composition evolution as a function of energy, for each size bin we take from the 
simulation the log{E) distributions of contributing mass groups weighted by the parameters pi and p 2 
and with weights Wk representing the relative efficiency to trigger the underground apparatus. The 
resulting distributions from different size bins are summed together, and so we eventually obtain the 
simulated energy spectra of the two basic components that reproduce the experimental data. These 
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Figure 1. Energy spectra for p and Fe (upper) Light and Heavy (lower) components, and their sum. 


resulting spectra for p-Fe and Light-Heavy components, further corrected for the detection efficiency of 
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the surface array, are reported in Fig. together with their sum. Within such procedure, for each bin in 
log(E) we calculate < logA > from the expression: 


< logA >= 


PiMlLog{AP) + p2WkM[^Log{A^‘^) 
Ml+p2WkMP- 


( 2 ) 


where and are the numbers of simulated events from p and Fe contributing to the k-th energy 
bin. The same expression, with obvious modifications, is used for the Light-Heavy composition. The 
results are depicted in Figs. the two curves delimiting the error band. The results show an increase in 
< logA > through the energy corresponding to the knee position. 


3. Conclusions 

The analysis of the Ng-N^ >1.3 TeV) data collected by the MACRO/EAS-TOP collaboration at 
the Gran Sasso Laboratories points to a primary composition becoming heavier around the knee of the 
primary spectrum (in the energy region 10^® — 10^® eV). The result is in very good agreement with the 
measurements of EAS-TOP alone at the surface, wich detects muons with energy > 1 GeV, using the 
same (QGSJET) interaction model 0. GeV and TeV muons are produced in different kinematic regions: 
in the central one and at the edges of the fragmentation region respectively. The present measurements 
show therefore a good consistency of the interaction model in describing the yield of secondaries over a 
wide rapidity region. 
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Figure 2. < logA > vs energy for p/Fe (upper) and Light/Heavy (lower) compositions. The two lines 
include the statistical uncertainities. 



















































